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Abstract. We study the apparent discrepancy between the standard PQCD predictions for the meson and 
baryon ratios and muhiplicities at moderate high pt > 2 GeV and recent experimental measurements in 
Au + Au collisions at \/sjvjv ~ ^30 GeV at the Relativistic Heavy Ion Collider (RHIC). We show that the 
differences, most pronounced in central collisions, can be explained by a strong non-perturbative baryon 
Junction component, which dominates the currently accessible experimental pt window and the non-abelian 
energy loss of fast partons propagating through hot and dense medium. The recently introduced two com- 
ponent hybrid model, which combines a quenched jet PQCD calculation in the Gyulassy-Levai-Vitev (GLV) 
formalism and a phenomenological "soft" part, is further elaborated to take into account the full 3D expansion 
in the pre-hadronization phase and include particle flavor dependent "soft" inverse slopes as suggested by the 
baryon Junction picture. We show that such approach can resolve what seems to be a factor of ~ 2 difference 
in the moderate high px suppression of tt" and h~ as recently reported by the PHENIX collaboration. The 
observed quenching of the high pt particle spectra and the large p/tv~ and p/iv^ ratios as a function of pT 
are found to be consistent with a creation of a deconfined phase and non-abelian energy loss of fast partons 
in a plasma of initial gluon rapidity density dN^ /dy ~ 1000 



INTRODUCTION 

The preliminary data from the first Au + Au run of the Relativistic Heavy Ion Collider (RHIC) at 
y/i^^ = 130 GeV was reported by the experimental collaborations (BRAHMS, PHENIX, PHOBOS and 
STAR) at the Quark Matter 2001 [1] conference. It suggests the observation of qualitatively new phenomena 
in comparison to the CERN Super Proton Synchrotron (SPS) Pb + Pb program results at \/sjvAf ~ ^'^■'^ GeV. 
First, the moderate highpx spectra of neutral pions in the 10% central collisions were reported by PHENIX [2] 
to be suppressed by a factor 2-4 relative to the PQCD predictions scaled by nuclear geometry {TabO^), the 
Glauber profile density per unit area at an impact parameter b). Second, the analysis of semi-peripheral 
collisions by STAR [3] shows that the elliptic fiow V2{pt) is compatible with the monotonically growing 
hydrodynamic predictions for small transverse momenta, but saturates at ~ 0.15 above px > 2 GeV. Third, 
studies of the charged particle rapidity density as a function of centrality by PHOBOS and PHENIX [4] 
revealed the expected (but previously not observed) TabO^) ^ ^pirt niini-jet component [5]. Finally, in the 
Pt = 1 — 3.5 GeV window where both 7r° and h~ have been measured, the 7r° are quenched by a factor of 2 
more than the negative hadrons. 

While the data are preliminary, the analysis of Ref. [6] suggests that these first results from RHIC may 
be interpreted as evidence for the production of a dense gluon plasma in central and semi-central reactions. 
The key diagnostic tool for probing the properties of the medium in that analysis is jet quenching [7]. 
Estimates based on high px suppression of negative hadrons indicate that the initial gluon density may 
have reached pg{T ^ 0.2 fm/c) ^ lOOpoj where po = 0.15 fm^"^ is the nuclear saturation density. At RHIC 



energies, \/sjvAr — ~ ^00 GeV, the initial flux of high pT jets becomes sufficiently large and the produced 
gluon plasma is sufficiently dense and long lived to make jet quenching finally observable. In contrast, at 
SPS energies the somi-hard tt" spectrum reported by the WA98 collaboration [8] was found to be actually 
enhanced by a factor of 2 at px ~ 2 — 4 GeV due to the Cronin effect [9]. We emphasize the importance of 
studying p + A systems to resolve the strength of the Cronin versus nuclear shadowing effects on the high 
Pt yields. 

An unexpected result at RHIC reported by PHENIX [10] was that the flavor composition of high pt > 2 
GeV negative and positive pt differential hadron yields may actually be dominated by anti-protons and 
protons respectively. Such behavior has not been observed so far in pp or AA collisions. In the conventional 
QCD multi-particle production phenomenology high hadrons arise from the fragmentation of colored 
partons, {q, g), whose spectra are computable via PQCD. If hot and dense matter is created in AA reactions, 
those partons naturally lose energy through elastic and radiative multi-particle final state interactions. The 
quenching of high pt partons is thus expected to result in the quenching of all their hadronic fragments 
including p and p. 

To resolve the new high pr baryon puzzle at RHIC we propose that the non-perturbative baryon Junctions 
component [11] of hadronic spectra provides a natural mechanism that could account for the anomalously 
large baryon/meson ratios. This component competes with the semi-hard PQCD processes even in pp and 
pA reactions but has only now been exposed in central Au + Au at RHIC, because the usual perturbative 
QCD pion and kaon component is strongly suppressed due to jet quenching. The Junction mechanism is 
further supported by the observation that in central collisions at midrapdity p/p « 0.65 and dN^ /dy ~ 
dN^/dy — dN^/dy ~ 15 [12], indicating a very significant baryon number transport (stopping) five units of 
rapidity from the fragmentation regions. 



HARD AND SOFT HADRONIC COMPONENTS 

The PQCD approach expresses the differential hadron cross section in p+p^h + X through a con- 
volution of the measured structure functions f{xa,Qa)a/p for the interacting partons (a = a,b), with the 
parameterized fragmentation function Dh/c{z,Qc) for f^e leading scattered parton c into a hadron of fla- 
vor h and the elementary parton-parton cross sections d(7'^°-^~*'^'^^ jdi. Next to leading order corrections are 
accounted for by the constant K ^ 1. Other degrees of freedom for PQCD fits to data include the choice 
of the renormalization scale and the parton broadening due to initial state radiation, which can be in- 
cluded through a normalized fcx-smearing distribution. A simple functional form is the Gaussian smearing 
/(/ct) = e~'^T/{fcT>/7r^fc2 where (fc^) — 0-6 — 1-0 GeV^ in pp, pp collisions. In AA reactions one expects 
larger initial parton broadening due to fcx kicks. In the absence of a medium, the invariant hadron inclusive 
cross section is given by 

^^~dfp ^ J dZcdXadXb J d'^'kTad'^k.Tbfi^Ta)fO<-Tb)fa/p{Xa,Ql)fb/p{Xb,Ql) 
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Dh/c{Zc,Qt) — -I 5{s + u + t), (1) 

where Xa,Xb are the initial momentum fractions carried by the interacting partons, Zc = Ph/Pc is the 
momentum fraction carried by the observed hadron. 

In the case of jet production inside a hot and dense medium, the fast parton loses a fraction of its energy 
due to gluon radiation induced by final state interactions. The full solution for the gluon radiative double 
differential distributions for the case of finite jet energies and few collisions was found in [7] in a form easily 
integrated numerically (to AE) taking into account both the kinematic bounds and the running coupling 
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FIGURE 1. In Fig. la the gluon jet fractional energy loss in the GLV formalism is shown for three 
initial gluon rapidity densities dN^ /dy — 200,400, 1000 for both ID and 3D expansion. The asymptotic 
analytic approximation ( shown here for dN^ /dy = 400) is inapplicable at RHIC. Fig. lb shows the ratio 
of charged particles computed in the soft+GLV quenched model relative to the conventional PQCD that 
accounts for the UAl data. Also shown is the greater quenching predicted for tt" with dN^ /dy = 400. 

We have also studied the effects of a strong transverse flow in addition to the boost-invariant longitudinal 
Bjorken expansion [13] and for the azimuthally averaged energy loss and found the simple result: 

^Ef2 ~ Y^^AE^^, z„ax - SR , (2) 

where v± is the average transverse expansion velocity and Zmax = is the duration of the energy loss phase 
{S < 1/(1 — v±)). While the azimuthal asymmetry of AE is sensitive to the details of the expansion, the 
mean is not. For wj_ = 0.6 and /3 = 1.5 AE^^ ~ 80%Ae[^. Gluon jet mean fractional energy loss for three 
typical plasma densities is shown in Fig. la. Comparison to the asymptotic result signifies the importance 
of finite kinematics. Non-abelian energy loss affects hadronization through rescaling the hadron momentum 
fractions Zc and modifying the fragmentation functions in Eq. (|^) [6]. We note that because of the GLV 
energy loss AE/E oc pg is very weakly energy dependent (see Fig. la), the effect of energy loss fluctuations, 
can be absorbed as a gluon density renormalization by a factor ~ 2 — 2.5 [14]. 

The PQCD fragmentation picture gives adequate description of particle multiplicities in pp and pp col- 
lisions for > 2 — 3 GeV but over-predicts soft hadrons due to its power law divergence. We use the the 
results of the soft string models that fit pp and pp data to parametrize by a simple exponential behavior 
(in the m limit) the low pt region. If similar soft string dynamics drives pp, pA and AA collisions 
then most of the global observables with the exception of the baryon-antibaryon asymmetry can also be 
inferred from already existing data. The FNAL E-735 data [15] on hadronic spectra in pp at ^/s = 200 GeV 
suggests values for the particle ratios that come very close to preliminary RHIC data (shown for comparison 
in parentheses). The p/n = 0.07 (0.08). For K/n = 0.11 (0.15) we find ~ 30% deviation. For pions and 
kaons one infers from charge and strangeness conservation that 7r+/7r~ = 1.0 (1.0) and / = 1.0 (1.1). 
The ratio that cannot be fixed from data on net baryon-free pp collisions is p/p. However, it can be predicted 
by baryon transport (Regge) theory, i.e. dN^ /dy — Z/?cosh((l — a_B(0))j/)/sinh((l — aB(0))yniax) with a 
coupling f3 = 0{1). In the baryon Junction picture q;b(0) ~ 1/2 [11]. For Au + Au collisions at RHIC at 
y/s = 130 AGeV (Ymax = 5.4), the prediction dN^/dy ~ dNP/dy - dNP/dy ~ 11 is close to the preliminary 
data. The predicted Junction trajectory slope, a'j ~ implies that the effective string tension k' is three 

times higher than k = l/(27ra'^) « 1 GeV/fm. This naturally leads to ~ 3(px)h, i.e. Tg ~ V^Tq . 

Given the pion mean inverse slope, Tq ~ 220 MeV, it follows that the proton and antiproton slopes should 
be Tq ~ 400 MeV for the junction (high pt) component of those yields. 




FIGURE 2. The p/tt^ and p/ir^ ratios are shown in Figs. 2a, 2b as a function of px for 3 different 
impact parameters b = 0, 6, 12 fm. The ratio of out fits to prehminary minimum bias PHENIX data is 
shown for comparison. 

Published RHIC data [4] constrain the charged particle multiplicity at midrapidity to dN^/dy{h = 0) ~ 
600. Thus with one parameter, the pion mean inverse slope (the — 275 MeV is irrelevant for the study of 
the P/tt^ and p/n^ ratios), the above theoretical analysis based on the string and baryon Junctions picture 
results in a solvable set of equations for the parameters in the phcnomcnological soft component of the 
hadronic spectrum: 



(3) 



This component is assumed to scale with the number of participants, although the Junction component may 
scale faster [11]. Eqs. (|lj|^) give the hadron inclusive multiplicity 

dyd'^PT dyd'^pT dyd^pT ' 

where TaaO^) is the Glauber profile density per unit area. The power law PQCD component is limited to 
the range px > 2 GcV. If the Junction and Junction loop components scales nonlinearly with Npart (b) , as in 
HIJING/i?i3, then the rate of change of the maximum oi Rb with increasing b would change correspondingly. 



CONCLUSIONS 

Fig. lb shows the ratio of the charged particle multiplicities computed in our model relative to unquenched 
scaled PQCD (which fits to the UAl data [16]). The high px tails show approximately constant suppression 
driven by the initial gluon density. Similar computation is shown for neutral pious for mean energy loss 
driven calculation with dN^ /dy — 400. In the intermediate = 2 — 5 GeV region the factor of 2 difference 
in the apparent suppression arises from the anomalously large baryon and antibaryon component predicted 
by the Junction picture. 

In going from central to peripheral collisions the smaller partonic densities result in larger typical sizes of 
strings and Junctions and somewhat smaller mean inverse slopes of the resulting soft hadrons. This effect is 
taken into account by introducing a difference AT^ = -10 MeV, AT^ = -25 MeV and AT^ = -50 MeV 
between central b = fm and peripheral b = 12 fm collisions. The major distinction between various 
centralities, however, is the different strength of jet quenching due to the different lengths of the path that 
the jet travels inside the plasma and the different density of the medium. Thus in peripheral collisions the 
baryon/meson ratio is expected to stay below unity similar to the pp case. Figs. 2a and 2b show the hybrid 
two component model predictions for the p/ir^ and p/n^ ratios as a function of transverse momentum and 



centrality. While the precise values Rb max of the maximum baryon and antibaryon excess are sensitive to 
the uncertainties in the estimates of the plasma density and the mean inverse slopes, the characteristic 
and b behavior is not. In Figs. 2a and 2b we have included the ratio of our exponential fits in a narrow 
Pt — I — 2.1 GcV window to the preliminary PHENIX minimum bias data [10]. The error bars reflect the 
one a uncertainties in the slope and intercept parameters. 

Unlike hydrodynamics the dual soft+hard hybrid model predicts that the baryon and antibaryon anomaly 
is in fact limited to finite p window and the baryon to meson/ratios eventually converge to the PQCD 
dominated base at high p^. It also predicts the systematic decrease of the predicted ratios with impact 
parameter which can be tested at RHIC. With higher statistics the dynamics of the baryon/meson ratios as 
a function of centrality may fact be probed in the current experimentally accessible range. 
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